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Abstract 

We propose a method for realizing true, real-space imaging of charge dynam- 
ics in a periodic system, with angstrom spatial resolution and attosecond 
time resolution. In this method, inelastic x-ray scattering (IXS) is carried 
out with a coherent, standing wave source, which provides the off-diagonal 
elements of the generalized dynamic structure factor, S(q 1 , q 2 ,u)), allowing 
complete reconstruction of the inhomogeneous response function of the sys- 
tem, x( x ii x 2, t). The quantity x has the physical meaning of a propagator 
for charge, so allows one to observe— in real time— the disturbance in the 
electron density created by a point source placed at a specified location, x± 
(on an atom vs. between atoms, for example). This method may be thought 
of as a generalization of x-ray crystallography that allows refinement of the 
excited states of a periodic system, rather than just its ground state. 
Keywords: 

1. Introduction 

Electron dynamics underlie all fundamental phenomena in chemistry, 



biology and materials physics. Recent advances in attosecond laser sources 
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have created widespread interest in studying such electronic processes in 
real time, particularly in cases where spatial information, e.g., about the 
detailed configuration of electron wave packets, can be inferred. [H El E] 

We recently proposed a different approach to attosecond imaging based 
on inelastic x-ray scattering (IXS).pE] Unlike more common, time-domain 
approaches that exploit state-of-the-art laser technology, in the IXS ap- 
proach an x-ray photon modulates the electron density of the system, and 
the ensuing response is measured the momentum and frequency domain. 
Spatial and temporal information is then achieved by solving the inverse 
scattering problem, which for typical experimental setups yields time res- 
olutions in the range of few attoseconds. |1] This method, which we will 
refer to here as "IXS imaging," should be considered a member of the same 
broad technique class as x-ray crystallography, which is routinely used to 
obtain real space images— for example of protein structures— with angstrom 
resolution. [5] The IXS approach to studying dynamics is somewhat more 
restrictive than laser-based approaches, in that it strictly allows one to ob- 
serve the time-evolution of the electron density of the system in response to 
a charged, point source. [1] However, it provides explicit, real-space images, 
and contains no intrinsic limit on the achievable time resolution: attosecond 
or even zeptosecond resolution is achievable. So far, IXS imaging has been 
used to study collective electron dynamics in liquid water, [6j [7] excitons in 
large-gap insulators, [8] and to measure the effective fine structure constant 
of graphene.[9] 

Up to now, however, there has been a limitation on the imaging as- 
pect of this approach. The objective in IXS imaging is to determine the 

charge response function (or "propagator" ) , x( x ii x 2,t), which physically 
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represents the amplitude that a point disturbance in the electron density at 
location x x will propagate to x 2 after elapsed time In reciprocal space 
this quantity is a function of two momenta, x{Qii Q21 u ) ■> but conventional 
IXS probes only its diagonal components, x(Qi ~ 9> The Fourier trans- 
form of the latter quantity, x( x jt)-> is causal and quantitatively accurate in 
its time evolution. However, it depends on only one spatial variable, and 
corresponds to the complete response averaged over all source locations. [10] 
If the system of interest is homogeneous, for example a free electron metal 
such as aluminum, this spatial averaging results in no loss of information. 
However, in an inhomogeneous system, important, local features may be 
averaged out. 

To overcome this limitation one must measure the off-diagonal momen- 
tum components of x{Qii Q21 u ) 1 where q x 7^ — q 2 . It was shown many 
years ago that this is possible, at least in principle, by using coherent stand- 
ing waves. [TT| IT2] . In this approach, an x-ray standing wave is created by 
exciting a Bragg reflection in a crystal (or by using an external coherent 
source, such as a Bonse-Hart interferometer). Within a coherence volume, 
the x-ray photon lies in a quantum superposition of two distinct momenta, 
ki and k 2 . One then places the IXS detector (typically a backscattering 
analyzer) at some scattering angle, which defines two momentum trans- 
fers, q x and q 2 . Under these conditions, due to interference between the 
two scattering channels, the cross section includes terms that depend on 

X(<?1,<?2> W )- 

Could this approach be used to solve the averaging problem in IXS 
imaging? It is not obvious that it can. In standing wave techniques one 

1 By reciprocity, the probability of propagating instead from X2 to Xi is the same. 
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does not have complete freedom to choose the two momenta, ki and k 2 . 
Because the wave field is created by a Bragg reflection, the momenta are 
always related by a reciprocal lattice vector, i.e., ki~k 2 = G. Hence, one's 
access to reciprocal space is highly constrained, and it is unclear whether 
enough information is accessible to permit a full refinement of x( x i, x 2,t). 

In this article we examine this issue by analyzing a simple model of 
a single quantum particle in a periodic potential. We find that whether 
enough reciprocal space is accessible depends on the dimensionality of the 
problem. We show that the standing wave approach fails in both one and 
two dimensions, because the dimensionality of the accessible momentum 
space is lower than the dimensionality of the set of data needed to perform 
a refinement. In three dimensions, however, the technique should in prin- 
ciple be viable: by performing IXS measurements under a sufficiently large 
set of standing wave conditions— each defined by a distinct reciprocal lat- 
tice vector— it should be possible to reconstruct the complete, unaveraged 
response x{ x \i x 2, t). In this sense, this approach can be thought of as a 
new type of x-ray crystallography that allows refinement of the collective 
excitations of a periodic system, rather than just its ground state. 

2. Background 

We begin by reviewing spatially-averaged IXS imaging, as it has been 

implemented previously [HI El |9] (for a more thorough review, see ref. jl]). 

In IXS, a monochromatic beam of x-rays is impinged on a system, from 

which it scatters in all directions. An energy-resolving detector is placed 

at some angle and measures the spontaneously Raman scattered photons. 

This detector defines a transferred energy, u, and transferred momentum, q, 
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and the intensity measured is proportional to the dynamic structure factor, 
S(q,u), which is the Fourier transform of the density-density correlation 
function [Bill]. 

While S contains a great deal of information about excited states, it is 
of rather limited interest from the point of view of dynamics. As an auto- 
correlation function, the dynamical information contained in S is indirect, 
reflecting only the general length and time scales over which excitations take 
place. Fortunately, S(q,u>) is related to the response function, x(<7i>92> a; )> 
which reveals the real space and time dynamics, by the quantum mechanical 
version of the fluctuation-dissipation theorem, 

S(q,u) = ~ l _ e -pn u Imxfa-g,")- (1) 

Hence, in principle, IXS provides access to the true, causal dynamics of the 
system. 

Of course, S does not supply the entire \ function, but only its imaginary 
part. This restriction is the IXS rendition of the well-known phase problem 
in x-ray crystallography. [5] In the present case, the phase problem may be 
solved by recognizing that x is a causal function, i.e., the real part can be 
obtained from the Kramers-Kronig (KK) relation 

where V denotes the principal part. Inverse Fourier transforming gives a 

function, x(a3,t), that corresponds to the complete response, x( x i^ x 2,t), 

averaged over all space with the difference x = X\ — x 2 held fixed. [TO] In 

a homogeneous system, x{ x i t) represents a complete parameterization of 

the response. Fig. [lj for example, shows an image of the electron density in 

graphite 400 attoseconds after the source. [9] 
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Figure 1: x{ x i t) &t t = 400 as determined from a spatially-averaged IXS study of single 
crystal graphite, reproduced from [5]. 

Note that Eq.[2]impl ies a need for a significant quantity of experimental 
data. The response x(Qi ~Qi w ) is a strong function of the three-dimensional 
momentum, q, and must be sampled with enough range and point density 
to define an inverse spatial Fourier transform. Moreover, at every value of 
q, a spectrum must be obtained over a sufficiently large range to carry out 
the KK transform, Eq. [2j The number of unknowns can often be reduced by 
examining specific projections of reduced dimensionality [H] or by exploiting 
crystal symmetry [9] (see section VI). But, in general, this method requirs 
a large amount of data, and has only become possible because of the advent 
of high-brightness, third-generation synchrotron x-ray sources. 

3. IXS with standing waves 

To obtain unaveraged images, i.e., that will give a faithful representation 
of the dynamics of systems that are inhomogeneous, the previous approach 
must be extended to determine the off-diagonal response, x(QxiQ2i cj )- As 
discussed in Section 1, this might be accomplished by the use of an x-ray 
standing wave. The most straight-forward way to produce a standing wave 
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Figure 2: (a) Layout of a standing wave IXS imaging experiment. Two identical crystals 



of the system of interest are mounted nondispersively. The first, which is asymmetrically 
cut, collimates the beam, and the second generates the standing wave. k\ 1 ki, and k% 
are the momenta of the incident, Bragg reflected, and inelastically scattered photons, 
respectively. ^ describes rotations of the second crystal around G. |(b)| Standing wave- 
field formed by incident photon fci and Bragg reflected photon Nodal planes of the 
standing wave-field are highlighted in blue. The nodal planes are shown to be coincident 
with the lattice planes, but the phase of the standing wave-field, and hence the position 
of the nodes, is tunable. 
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is by exciting a Bragg reflection. Hence, we will from this point assume that 
the system of interest is available in the form of high quality single crystals, 
such as those commonly used in structural biology. [5] 

The proposed geometry for the experiment is illustrated in Fig. [2j Two, 
identical crystals of the system of interest are mounted in the geometry 
shown, with the first cut asymmetrically with respect to a particular set of 
lattice planes, defined by reciprocal lattice vector, G, and the second cut 
symmetrically. As described by dynamical diffraction theory, [H] the first 
crystal compresses the angular divergence of the beam so that it is much 
less than the Darwin width, 9\y, of the second, symmetric reflection. Under 
these circumstances, the second crystal will generate a coherent standing 
wave, in which the photon lies in a quantum superposition of two momenta, 
described by the quantum state [UJ H2J US] 

\^ = (9A iai +92a{ 2a2 e^)\m). (3) 

In this expression gi and g 2 are the amplitudes of the incident and diffracted 
beams, respectively, and 7 is the phase shift between the two. cika anni- 
hilates a photon with momentum k and polarization state a, and \m) is 
the initial many-body state of the electron system. Note that, because the 
wave field is created by a Bragg reflection, k\ — k 2 = G. Nevertheless, be- 
cause Bragg scattering is an elastic process, u(ki) = uj{k 2 ). For notational 
simplicity, we take lo\ = u(kx) = u(k 2 ) and u 3 = u(k 3 ) (Fig. 2) . The 
parameters g 2 and 7 may be controlled by fine adjustment of the angle of 
either crystal, and both quantities may be accurately calculated from dy- 
namical diffraction theory. For a thorough review of x-ray standing wave 

techniques, we refer the reader to Refs. [HJ ITS] . 
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Once the standing wave is established, the initial state in the inelastic 
scattering process is that shown in Eq. |3j The final state is an ordinary 
plane wave state, given by 

l/> = *L»- ( 4 ) 

where k 3 is the wave vector of the inelastically scattered photon, and \n) is 
the final many body state to which the electron system was excited. 

The scattering from \i) to |/) is mediated by the usual interaction be- 
tween light and matter. The Hamiltonian for a system of non-relativistic 
electrons interacting with the electromagnetic field is [T7] 

H = H Q + ^- J i) ] A-pi)dx + ^-^ J pA 2 dx (5) 

where Hq is the many-body Hamiltonian of the electron system, ip(x) is 
the electron field operator that annihilates an electron at position x, and 
p(x) = $ (x)ip(x) is the electron density operator. A is the quantized 
vector potential, 

A(x) = (^j 2 J2 ^= (tiA a e~ lk x + e ka a ka e lkx ) , (6) 

which either creates or annihilates a photon at location x. 

To lowest order in perturbation theory, only the second interaction, 
flint — {e 2 /2mc 2 ) J pA 2 dx, leads to scattering. [T3] The first order scattering 
amplitude is given by 

(f\ flint K) = — 7===(es ■ ei) (n\ p(-q x ) \m) 

2irhe 2 g 2 e^ A 
mV ^oo 3 u 2 
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where we have defined the two momentum transfers q 1 = k\ — k 3 and 
q 2 = k 2 — fe 3 . The transition rate IV^ at finite temperature T is given by 
Fermi's golden rule 

iw ~ (^) E 6 - W^i ■ ^\ 2 \ H p(«i) i™> i 2 

^ ' n,m 

+92^1 ■ ^2| 2 | (n\p(q 2 ) \m) \ 2 

+9i92e^(e 3 ■ e*)(eg • e 2 ) (m| p^) |n) (n| p(-g 2 ) l m ) 
+^i^e~ i7 (e* • ei)(e 3 • e 2 ) (m\ p(q 2 ) \n) (n\ p{-q x ) \m)] 
x5(E-E n + E m ), (8) 

where we have used the fact that the electron density is real, i.e., p\q) = 
p(—q). Here b m = e~^ Em /Z is the Boltzmann factor and E is the energy 
transferred to the sample by the photon. 

The quantity that is relevant to the experiment is the doubly-differential 
scattering cross section [TT] 

d 2 a 1 d 2 N 



dttdE $ -'OttOE- 

where d 2 N / OVLdE is the density of final states and $ is the incident flux. 
Because the final state contains a single photon with polarization a, we 
have d 2 N/dVtdE = Vu 2 /8n 3 hc s and $ = c/V. Substituting into Eq. M we 
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have 



fj_ \ — [^|63-ei| 2 | <ti| p( qi)\m} 12 

\ mc z J uji z — ' 

N ' n,m 



dE'dn 



+02 1 35 ■ g a| 2 | (n\p(q 2 ) \m) | 2 

+9i92e tl (e 3 ■ ej)(eg • e 2 ) (m\ p{q l ) \n) (n| p(-q 2 ) \m) 
+gig 2 e~ n (el ■ ei)(e 3 ■ e^) (m| p(g 2 ) |n) (n| p(-qfi) |m) 
x 5(E-E n + E m ). (10) 

By inspection, we recognize the first two terms as the usual dynamic struc- 
ture factor. [TBI 



S(q, lo) = J2 b m \ (n\ p(q) \m) \ 2 5(E - E n + E m ), (11) 

n,m 

which is the correlation function measured in normal IXS experiments. 



Hence, Eq. [TO] can be rewritten 

d 2 a f e 2 \ u 3 r 2 A 2 2 A 2 

+5 , i5 , 2e' 7 (e 3 ■ e*)^ ■ e 2 )S{q x ,q 2 ,u) 

+g 1 g 2 e^(et-e 1 )(e 3 -e* 2 )S(q 2 ,q 1 ,u J )] , (12) 

where u = u± — u 3 is the energy transferred to the sample, and we have 
defined a nondiagonal, generalized dynamic structure factor 

S(qi,q 2 ,^) = ^2b m (m\ p(q 1 ) \n) (n\ p{-q 2 ) \m) 5(hu - E n + E m ). (13) 

n,m 

This quantity physically represents the degree to which charge fluctuations 

with wave vector q l are correlated with those with wave vector q 2 . In the 

next section we will see that this generalized correlation function is related 

to the desired off-diagonal components of the susceptibility, x(<h) ~q2i 0J )- 
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4. Generalized Fluctuation Dissipation Theorem 

To proceed further we need a fluctuation-dissipation theorem that ap- 
plies to the generalized function, S(q 1 ,q 2 ,cj). The retarded response func- 
tion, x, is defined in real space as a quantum mechanical thermal average, [18 



X(xx, x 2 , t x - t 2 ) = ~ ^2 b ™ H \p{ x ^ *0> P( x *> l m ) " t2 )' ( 14 ) 

m 

where 6(t) is a step function. In reciprocal space the imaginary part is given 
by 

Imx(Pi,p 2 ,w) 

= - 7r^6 m [(m\ p{p x ) \n) (n| p(p 2 ) |m) <5(fru; - £ n + E m ) 

n,m 

- (m| p(p 2 ) |n) (n| p(p x ) |m) <5(ftu; + £ n - £ m )] . (15) 
where the \m) are, again, many-body eigenstates of the electronic system. 



Eqs. 15 and 13 together imply that 



Imx(<?i,<? 2 > w ) 
= - nS(q X i -q 2 ,w) 

+ 7r^6 m (m|p(q 2 ) Hp(<?i) \m) 5(-huj - E n + £ m ). (16) 

Switching the indices n and m in the second term and noting that b n = 
bme' 13 ^, we arrive at 

This expression is a generalized form of the fluctuation-dissipation theorem 

that applies to a system that lacks translational symmetry. This result is 
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significant because it implies that the cross section Eq. [T2j which contains 
terms involving the generalized S, provides access to the needed off-diagonal 
terms of \- Of course, these terms are tangled up with the usual diagonal 
response, but they may be isolated by varying the phase shift, 7, subtracting 
the diagonal structure factor terms, and properly normalizing. [Til IT5] 



While the general form of Eq. 17 is not unexpected it reveals a crucially 



important detail, which is that there is a sign difference on the second 
momentum, q 2 , between x an d S. As we will see below, this sign is crucial 
for the viability of standing wave IXS as a technique. 

5. Quantum particle in a periodic potential 

We now address the question of whether standing wave IXS can pro- 
vide, at least in principle, enough information to carry out a refinement of 
the full real space response, x{ x ii x 2, t). We begin by pointing out that, 



because of the negative sign on q 2 in Eq. 17, the experimental standing 
wave constraint, k\ — k2 = G, is not as severe as it may seem. Because the 
system is periodic, there is also a constraint on the response itself, i.e., 

X(x ll x 2l t) = x(xi + R,x 2 + R,t), (18) 

where R is a real space Bravais lattice vector. In reciprocal space this 
constraint has the form j 10 



x(Pi,P2>u)=x(Pi,G-Pi,w) ( 19 ) 
i.e., x is nonzero only when p 1 + p 2 = G. In terms of the actual momenta 



in the experiment, related by Eq. 17, this constraint implies that G 
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Pi + P2 = Qi ~ Q2 = k\ — — k2 + k% or, in other words, that x is nonzero 
only when 

k 1 -k 2 = G. (20) 

This is, however, precisely the same as the constraint on the momenta from 
the standing wave condition itself. 



In other words, the minus sign on q 2 in Eq. 17 has the consequence that 
the standing wave constraint is no constraint at all. While the standing wave 
approach does not allow access to all possible combinations of momenta, q x 
and q 2 , not all combinations are needed; the response function, x is nonzero 
only for a select subset of momenta, and this subset is exactly what is 
accessible in the experiment. 

The last issue is whether, at a fixed value of G, there is sufficient flexi- 



bility in the remaining momentum (p 1 in Eq. 19) to perform a refinement 
of the complete x- To answer this question requires a specific, microscopic 
model. For this purpose, we consider the case of a single, spinless quantum 
particle traveling in a periodic array of harmonic wells, 

V(x) = -vr\ + v \ x ~ R] \ d ( r o ~\x-R\), (21) 

R 

where r is the radius and v the depth of each well. For simplicity we assume 
the system is simple cubic with lattice parameter a, and that r < a. While 
this model is not particularly physical, the conclusions we draw from it 
about momentum constraints will be quite general. 

If the wells are deep, the lowest energy band can be computed in the 
LCAO or "tight binding" approximation. Its dispersion is 

oj(k) = ujq — 2th [cos(k x a) + cos(k y a) + cos(k z a)] — vr\ (22) 
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where ooq is the ground state energy of an isolated well and —th is the 
hopping parameter between wells. The corresponding bloch waves have the 
form ij)k{x) = J2 R 4>(x — R)e' tk ' R /\^N, where the tight-bonding orbitals are 
just the ground state of the simple harmonic oscillator, 

<j){x) = (na 2 yh~^ (23) 

and N is the number of wells. The electron field operator is then given by 

4>(x) = ^c k ^ k {x) (24) 

k 

where c k annihilates an electron in state ip k , and d is the dimensionality 
(d=l, 2, or 3). 

We now wish to compute the generalized dynamic structure factor, 
S(q 1 , q 2 , ou), for this model, and determine— for any given value of G— whether 
it is experimentally possible to access a sufficient range of momenta. Be- 
cause the bands disperse, the scattered intensity depends not just on the 
momenta, but also on the value of the transferred energy, u, resulting in a 
problem that is effectively seven dimensional. To simplify the discussion, 
we define an auxiliary quantity 

/oo 
du S(q 1 ,q 2 ,uj) 
-oo 

= h m H p(Qi) \ n ) ( n \ P(-<? 2 ) \ m ) , (25) 

n,m 

which describes the total amount of spectral weight to be found at a given 
combination of q 1 and q 2 , irrespective of the value of u. For simplicity, we 
take the zero temperature limit, bo — 1 and 6 m ^ =0, in which case 



where we have now labeled the states in terms of the momentum, i.e., \k) 
denotes the electron in eigenstate ipk( x ), |0) being the ground state. 

In momentum space, the density operator p(x) = fo(x)ip(x) has the 
form 

p(q) = H<i)52<ic k+q+G (27) 



k,G 



where $(<?) is the Fourier transform of |0(:e)| 2 . Substituting Eq. 
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into 



Eq. [26] gives 

£(«!, q 2 ) = NQfafifa) ^ 6 Ql ^ G (28) 

G 

where we have used the fact that |0(;c)| 2 is centrosymmetric, i.e., $(q) = 
$(— q). As expected, the intensity scales with the number of scatterers, 
N, and is nonzero only when q l — q 2 = G, which required by periodic- 
ity. Subject to this constraint, the intensity is just given by the quantity 
$(q 1 )$(q 2 )- We now evaluate where this product is nonzero and whether 
it can be adequately sampled in a physical experiment. 

5.1. One- Dimensional Case 

For illustration, it is useful to consider the problem first in one dimen- 
sion, before moving on to the physical cases of two and then three dimen- 
sions. 



In one-dimension Eq. [21] describes a chain of harmonic wells. This case is 
not experimentally realizable, since it requires the two momentum transfers 
to be collinear, which is impossible except in the limit of infinite energy. 
We can, however, draw several important conclusions about the functional 
form of £ from this case. 
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In one dimension the momenta q\ and q2 are scalars, and the reciprocal 
lattice vectors G = 2irh/a are indexed by a single integer, h. The resulting 
value of £(<7i, 52) is illustrated in Fig. [3j Because of the constraint q\ — q 2 = 



G, £ is nonzero only along the set of lines depicted in Fig. 3(a) 



Fig. 3(b) shows the magnitude of i(qi,q\ — G) plotted against q± along 
a selection of these lines. Two important observations can be made from 
this plot. First, the intensity along a given line is not uniform, but has the 
shape of a Gaussian centered half-way between qi = and q\ = G, whose 
momentum width Aqi = er -1 . Second, the intensity varies from one line 
to the next, decaying with increasing \G\, again with a Gaussian envelope 
with width a -1 . 

From this simple case we see an important relationship between the 
degree of inhomogeneity of the system and the range of momenta that must 
be sampled experimentally to reconstruct x{ x ii x 2, t). The narrower the 
wells, i.e., the more inhomogeneous the system, the larger the number of G 
values that must be sampled, and the larger the range of q\ that must be 
measured for each G. Hence, the number of G values that must be sampled 
is not infinite, but is of order a/ a, which is a measure of the strength of 
local field effects. For the parameters chosen here, one need only sample 
up to h = 2 to acquire > 99% of the spectral weight that is available. In 
the limiting case of a homogeneous system, a — >■ 00, £ is nonzero only for 
G = 0, and conventional IXS can completely parameterize the response, as 
expected. 
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Figure 4: Relative magnitude and orientation for the experimental geometry in the two- 
dimensional, co-planar case. Here G = 2tt(2, 0)/a. 
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5.2. Two-Dimensional Case 

In two dimensions, a standing wave experiment is physically realizable, 
and corresponds to all the rays in Fig. [2] lying in a single plane. In this case 



Eq. 21 describes a square, planar array of harmonic wells, the quantities 
q l and q 2 are two-component vectors, and the reciprocal lattice vectors 
G = 2ir(h, k)/a are described by two indices, h and k. 

The quantity q 2 ) is now a function in four- dimensional space and is 
not easily illustrated. Extending reasoning from the one-dimensional case, 
however, we expect £ to be nonzero only along discrete, two-dimensional 
sections through this space, again defined by the constraint q x — q 2 = G^\ 
Moreover, we expect the magnitude of £ to decrease with increasing \G\, 



and to be substantial only for y/h 2 + k 2 < 2. By analogy with Fig. 3(b) 



we plot in Fig. [5] several sections of £(<Zi,<Zi — G) for selected values of G. 



As in Fig. 3(b) £ is a Gaussian function centered at the midpoint between 
q 1 = and q 1 = G, whose width is isotropic and equal to a -1 . 

Because the two-dimensional case is experimentally realizable, it is now 
possible to address the fundamental question (see Section [T]) of whether 
enough information is accessible with standing wave IXS to reconstruct 
x{xi,X2,t). To do so, one must experimentally parameterize a set of two- 
dimensional surfaces, such as those illustrated in Fig. [5j each defined by a 
distinct G. 

For a given choice of G, the angle of the beam with respect to the 
crystal is fixed by the Bragg condition. Hence, the only means of adjusting 



2 Generally speaking, in d dimensions £ is a scalar function residing in a 2d-dimensional 
space, and the constraint q 1 q 2 = G defines a discrete set of d-dimensional sections 
through this space. 
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the momenta is to adjust the scattering angle of the IXS analyzer. Doing so 
has the effect of tracing out a circular trajectory in momentum space, that 
has radius k% and intersects both the origin and q 1 = G. For illustrative 
purposes we display in Fig. [5] (black lines) each of these trajectories, for the 
specific case of an x-ray beam energy E = 15.5 he/ a. 

A serious problem is now evident. One must experimentally parameter- 
ize two-dimensional surfaces such as those shown in Fig. [5j however the set 
of experimentally accessible points is only one- dimensional, corresponding 
to a cut through each of the needed surfaces. Hence, we are faced with a 
crisis of dimensionality: The set of accessible information is of lower dimen- 
sionality than what is needed to refine the response function. 

The dimensionality could be increased, of course, by adjusting also the 
beam energy, as illustrated in Fig. [6] This changes the radius of the circular 
section, in principle allowing one to sweep out a two-dimensional surface. 
Unfortunately, the nature of IXS is such that each incident energy (apart 
from modest amounts of tunability) requires a distinct experimental setup 
with a distinct energy analyzer. We conclude that the two-dimensional case 
of standing wave IXS imaging, while experimentally realizable, lacks the 
momentum flexibility to be viable for imaging. 

5.3. Three-Dimensional Case 



Finally, we address the three-dimensional case, in which Eq. 21 describes 
a cubic lattice of harmonic wells. All momenta now have three components, 
and the G vectors are described by three integers, h, k, and I. £(<?i,<7 2 ) 
is a function in six-dimensional space, but periodicity dictates that it is 
nonzero only along discrete, three-dimensional contours parameterized by 
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(b) G=% (2,0). 

-5 5 -5 5 

91, (A"') quiK 1 ) 

(c) G=f (-2,0) (d) G=f (0,-2) 

Figure 5: £(<7 l5 <?! — G) plotted for the two-dimensional case against the two components 
qi x and qi y for several values of G. £ has the form of a Gaussian centered at the half-way 
point between q 1 = and q 1 = G. The solid lines are the contours traced out by rotating 
the IXS analyzer through 360°. 



(a) G=^(0,2) 
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Figure 6: Contour shown in Fig. 5(b) for several values of the beam energy. Red, 
dashed green, and dotted blue contours correspond to E = 12.3 he/ a, E = 15.5 he/ a, 
and E = 20 A he/ a, respectively. Continuously tuning the beam energy in this manner 
can, in principle, trace out a two-dimensional surface, though doing so is prohibitively 
difficult experimentally. 

Qi ~ Q2 = G. By analogy with the previous two cases, we expect the 
function £(<7i,qfi — G) to be a Gaussian with width o~ x centered on the 
point q 1 = G/2. 

In terms of the measurement itself, two types of motions are now possible 
that were not available in either of the previous cases: (1) The analyzer may 
now be rotated in two directions, i.e., both parallel and perpendicular to 
the Bragg plane. (2) The sample may be rotated around an axis parallel to 
G, changing the sample angles while maintaining the Bragg condition. The 
latter motion is usually referred to rotation". 

The analyzer motion, which is now two-dimensional, traces out a spher- 
ical shell in q l space. As in the two-dimensional case, this sphere has radius 
k± and intersects both the origin and the point q 1 = G. As before, this 
section is of lower dimensionality than the space spanned by the function 
£(9i>9i — G). An additional degree of freedom is needed to make the di- 
mensionality of the measurement match the that of the £ function. 
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This degree of freedom is supplied by the \I/ rotation. Rotating the 
crystal around the Bragg vector (i.e., maintaining the Bragg condition) has 
the effect of sweeping the spherical shell around an axis defined by the line 
connecting the origin and q 1 = G. In this manner, the shell sweeps out 
a three-dimensional volume in momentum space. The volume swept out is 
a torus, with major radius \fk i — G 2 and minor radius k± 7 as illustrated 
in Fig. 7. A plot of the intersection of this torus with a constant-contour 
plot of £(<Zij<Zi — G) is shown in Fig. |8} We conclude that the crisis of 
dimensionality encountered in the two-dimensional case does not take place 
here, and that standing wave IXS imaging— in principle— should be a viable 
technique in the real world of three dimensions. 

While the dimensionality of the experiment is adequate, the three di- 
mensional case has some blind spots. As one might expect, information 
outside the torus shown in Fig. [8] is not accessible. This is the standing- 
wave manifestation of the well-known diffraction limit, which says that it is 
not possible in a scattering experiment to extract information about features 
smaller than half a wavelength. 

Second, and less intuitively, there is a blind spot in the center of the 
torus (see Fig. [7]) whose shape is a three-dimensional vesica piscis with axis 
length \G\ and major radius k\ — a/ ki 2 — G 2 . It remains to be determined 
whether this inner blind spot poses a serious limitation on using standing 
wave IXS for imaging. We note, however, that the radius of the blind 
spot goes to zero as k\ — > 00, so it can always be made arbitrarily small 
experimentally by working at sufficiently high beam energy. 
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Figure 7: Toroidal volume swept out by combined rotation of the two detector motions 
and the sample azimuth, ty, in the three-dimensional case. The center contains a blind 
spot whose shape is a vesica piscis swept through its major symmetry axis. This blind 
spot may be made arbitrarily small by increasing the beam energy. 




Figure 8: Intersection of the toroidal measurement volume shown in Fig. 7 with a 
constant-contour plot of £(qi, qi — G). 
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6. Conclusions 

To summarize, we have found that standing wave measurements are, 
in principle, a viable approach to overcoming the translational averaging 
problem [10] in IXS imaging. j6j E] If successfully implemented, this approach 
would reveal the complete density response, x{ x ii x 2,t), which describes the 
electron disturbance created by a source placed at any arbitrary location, 
in a spatially inhomogeneous but periodic system— typically with Angstrom 
spatial and attosecond time resolution. This technique can be thought of 
as a generalization of x-ray crystallography that allows refinement of the 
excited states of a periodic system, rather than just its ground state, and 
represents the maximum that can be learned by interaction of light with 
matter in the regime of linear response. This technique would be most useful 
for imaging excitations in very inhomogeneous systems, such as molecular 
crystals, in which local field effects are significant and transverse collective 
excitations, such as transverse plasmons, can be important. 

Analyzing a simple model of a single quantum particle in a periodic 
potential, we have shown that standing wave IXS imaging is an innately 
three-dimensional measurement, in the sense that both out-of-plane ana- 
lyzer motions and sample ^-rotations are required to achieve a complete 
data set. Such an experiment would require (at the minimum) two copies 
of the crystal of interest: one that is asymmetrically cut to collimate the 
beam, and a second to create the standing wave field and function as the 
"sample" . The scattering experiment would then require six rotation axes: 
a conventional three-axis sample goniometer, a single-axis, two-theta ro- 
tation supporting an energy-integrating detector for measuring the Bragg 

diffracted beam, and a second two-theta rotation, with both in- and out-of- 
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plane degrees of freedom, supporting a backscattering, IXS analyzer. 

One might think that the data collection time required to refine, for 
example, the three-dimensional x( x ii x 2,t) for a molecular crystal would 
be enormous, but this is not necessarily so. Specifically, suppose one were 
studying a molecular crystal whose ground state density had been refined 
with conventional x-ray crystallography, to a resolution providing N carte- 
sian voxels in each unit cell. Such a refinement would have required of 
order N independent Bragg measurements (unless other information about 
the structure, such as its symmetry, were known). To refine x{ x ii x 2, t) for 
this structure with the same resolution, because the response is a function 
of both source and observation coordinates, would require of order N 2 mea- 
surements at a single time slice. To construct the complete dynamics, then, 
requires N 2 N t measurements, where N t is the number of points in the time 
series. This number is, without a doubt, impracticably large in nearly all 
conceivable cases. 

The time considerations, however, need not be so severe. For example, 

if one fixes X\, i.e., if one decides before the measurement where to "strike" 

the molecule, the number of unknowns reduces to NN t . Further, as in 

x-ray crystallography, knowledge of symmetry further reduces the number 

of unknowns. For example, all four of the independent measurements in 

Fig. [6j which are related by symmetry, would reveal the same experimental 

result; this could have been anticipated ahead of time from the crystal 

symmetry. Other tricks, such as taking one- or two-dimensional projections, 

compromising the resolution by probing only small order G values, etc., are 

always possible in specific cases. Hence, we expect the data collection time 

to be no more intensive than the spatially-averaged IXS imaging, which has 
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already been demonstrated. 0, EJ [9] 
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